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Abstract 

One of the most important questions in studying quantum computation is: whether a quantum computer 
can solve NP-complete problems more efficiently than a classical computer? In 2000, Farhi, et al. (Science, 
292(55 16):472^476, 2001) proposed the adiabatic quantum optimization (AQO), a paradigm that directly 
attacks NP-hard optimization problems. How powerful is AQO? Early on, van Dam and Vazirani claimed 
that AQO failed (i.e. would take exponential time) for a family of 3SAT instances they constructed. More 
recently, Altshuler, et al. (Proc Natl Acad Sci USA, 107(28): 12446-12450, 2010) claimed that AQO failed 
also for random instances of the NP-complete Exact Cover problem. In this paper, we make clear that all these 
negative results are only for a specific AQO algorithm. We do so by demonstrating different AQO algorithms 
for the same problem for which their arguments no longer hold. Whether AQO fails or succeeds for solving 
the NP-complete problems (either the worst case or the average case) requires further investigation. Our AQO 
algorithms for Exact Cover and 3SAT are based on the polynomial reductions to the NP-complete Maximum- 
weight Independent Set (MIS) problem. 



1 Introduction 

A quantum computer promises extraordinary power over a classical computer, as demonstrated by Shor HI 
in 1994 with the polynomial quantum algorithm for solving the factoring problem, for which the best known 
classical algorithms are exponential. Just how much more powerful are quantum computers? In particular, can 
a quantum computer solve NP-complete problems [2l more efficiently than a classical computer? NP-complete 
problems are the "hardest" problems in NP — in the sense that if one can solve one NP-complete problem 
efficiently (i.e., in polynomial time) then one can solve all the problems in NP in polynomial time. In 2000, 
Farhi et al. El 01 proposed an adiabatic quantum algorithm as an alternative quantum paradigm to directly solve 
NP-hard optimization problems, which are polynomially equivalent to their corresponding NP-complete decision 
problems. Apparently, the same idea to the adiabatic quantum optimization, under a different name of quantum 
annealing, was first put forward by Apolloni et al. in 1988, see 0[6l and references therein for a history of the 
field. 
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1.1 Adiabatic Quantum Algorithm 

An adiabatic quantum algorithm is described by a time-dependent system Hamiltonian 

Hit) = (1 - s(t))H init + s(t)ft prob | em (1) 

for t G [0, T], s(0) = 0, s(T) = 1. There are three components of %{.); (1) initial Hamiltonian: %(0) = H- m \ t ; 
(2) problem Hamiltonian: 1~L{T) = % pro biem; and (3) evolution path: s : [0, T] — > [0, 1], e.g., s(t) = T-L(t) is 
an adiabatic algorithm for a problem if we encode the problem into the problem Hamiltonian H pro biem sucn that 
the ground state of % pro biem corresponds to the answer to the problem. The initial Hamiltonian H m \t is chosen 
to be non-commutative with % pro biem an d its ground state must be known and experimentally constructable, e.g., 
%initai = — YlieV(G) ^i a f- Here T is the running time of the algorithm. According to the adiabatic theorem, 
if T-L(t) evolves "slowly" enough, or equivalently, if T is "large" enough, which scales polynomially with the 
inverse of the minimum spectral gap g mm (the difference between the two lowest energy levels) of the system 
Hamiltonian, the system remains in the ground state of %(£), and consequently, ground state of%{T) = % pro biem 
gives the solution to the problem. 

This computational model is referred as the Adiabatic Quantum Computation (AQC). It has been shownQ 
8 ] that AQC is polynomially equivalent to conventional quantum computation (quantum circuit model). For 
the optimization problem, the problem Hamiltonian can be expressed as a diagonal matrix in the computa- 
tional basis. That is, let / pro biem : {0, l} n — > M be a cost function of the optimization problem such that 
the minimum of the / pro biem corresponds to the solution of the optimization problem, then the corresponding 
problem Hamiltonian % pro biem is the Hamiltonian with / pro biem as the energy function: namely, % pro biem = 
Sxe{o l} n /problem (x) \x) (x\ (which needs to be expressible in polynomial resources, such as Eq. (4) ). Hereafter 
we use the cost function and the energy function of the problem Hamiltonian interchangeably. It is worthwhile to 
emphasize here that given a problem, there can be many possible cost functions and thus many possible problem 
Hamiltonians for the same problem. This restricted model is referred as the Adiabatic Quantum Optimization 
(AQO) (which is no longer polynomially equivalent to the quantum circuit model). We remark that this distinc- 
tion between AQC and AQO was not made by Altshuler et al. (9l (They referred both as AQO). In this paper, the 
focus is on this restricted model. 

1.2 What Does "AQO Fails" Mean? 

The NP-complete problems that were initially proposed for AQO by Farhi et al.|]3]|D were 3SAT and a special 
case of 3SAT — Exact Cover 3 (EC3). As an example, they proposed a clause-violation cost function as the en- 
ergy function of the problem Hamiltonian. Namely, / pro biem {x) = number of clauses violated by the assignment x. 
This cost function (with perhaps constant difference) has been adopted by almost all the other adiabatic quantum 
computation works. In particular, van Dam and Vazirani ifTOl claimed that AQO failed to solve a family of 3SAT 
instances by showing that the AQO algorithm with this specific clause-violation cost function had the exponen- 
tially small minimum spectral gap. See Discussion for more discussion on some other similar claims([ 1 1 , 12]). 
Recently, Altshuler et al. [9) claimed that AQO failed for the average case of NP-complete problems by arguing 
this specific clause- violation cost function based AQO algorithm failed for random instances of EC3, because of 
the Anderson localization phenomenon. While the cost function proposed is a "natural" one, nevertheless, it is 
not the only possible one. Recall that, according to the formulation of an AQO algorithm, the requirement of 
the problem Hamiltonian is that the ground state corresponds to the solution. There can be other cost functions 
with the same minimum (solution), e.g., see the reduction below. That is, there are other problem Hamiltonians 
that have the same ground state but a different energy spectrum. Given a problem, there are three components 
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(initial Hamiltonian, problem Hamiltonian, and evolution path) that specify an AQO algorithm for the problem. 
A change in one component (e.g. problem Hamiltonian) will result in a different AQO algorithm for the same 
problem. So how to prove that AQO fails for a problem? When is it sufficient to argue for one specific cost 
function (and thus specific problem Hamiltonian) and generalize to all other possible problem Hamiltonians? 

1.3 NP-Complete Reduction Based Problem Hamiltonian 

In this paper, we make clear that the arguments in both van Dam and Vazirani [10] and Altshuler et al. (9j for 
their specific AQO algorithm do not generalize. We do so by concretely describing another problem Hamiltonian 
for the same problem to which their arguments no longer apply. The counter-arguments we provide are simple 
enough that researchers who are not familiar with AQC and/or do not fully understand the arguments in (9l [lOl 
can easily follow. Our problem Hamiltonians are based on the polynomial reductions - NP-complete problems, 
by definition, can be polynomially reducible to each other - to the Maximum-weight Independent Set problem 
(see e.g. ||2l), which is one of the well studied NP-complete problems. Perhaps, a more interesting and important 
question is: would the NP-complete reduction make a difference for the adiabatic complexity of the problem? 
Recall that the reduction requires only the solution to be preserved, i.e. there is a polynomial time algorithm that 
maps the solution to the reduced problem to the solution to the original problem and vice versa (see e.g. |fT3l ). In 
other words, the reduction might only preserve the solution (i.e. the ground state) and alter the energy levels of 
the problem Hamiltonian. In Ifl4l . we demonstrate that even in small examples, the minimum spectral gap can be 
increased drastically when the excited energy levels are changed, due to the freedom in selecting parameters in 
the problem Hamiltonian. Thus, different reductions are possible, giving rise to different problem Hamiltonians, 
and thus different AQO algorithms, for the same problem. What are the time complexities of these different AQO 
algorithms? Unfortunately, at this point, the analytical analysis of these algorithms, which requires bounding the 
minimum spectral gap of the system Hamiltonian, remains challenging. Whether the time complexity of each of 
these AQO algorithms is polynomial or exponential remains open and requires further investigation. 

2 Methods and Results 

In the following, we first recall the NP-complete Maximum-weight Independent Set (MIS) problem, and a prob- 
lem Hamiltonian for solving MIS. Then we describe a simple reduction from Exact Cover to MIS that results in a 
different problem Hamiltonian (and thus different AQO algorithm) for Exact Cover. We then make clear that the 
argument in Altshuler et al. [9 ] does not apply to this algorithm. Similarly, we recall the well-known reduction 
from 3SAT to MIS, and point out that the argument in [ 10] does not apply to the MIS-based AQO algorithm. 

2.1 Maximum- weight Independent Set 

Formally, the Maximum-Weight Independent Set problem is as follows: 

Input: An undirected graph G(= (V(G), E(G))), where each vertex i G V(G) = {1, . . . , n} is weighted by 
a positive rational number q 

Output: A subset S C V(G) such that S is independent (i.e., for each i, j G V(G), i / j, ij G" E(G)) and 
the total weight of S (= YlieS c *) * s maximized. Denote the optimal set by mis(G). 

This is referred as the optimization version of MIS. The corresponding decision problem is: 

Question: Is there an independent subset S C V(G) such that the total weight is at least k? (where the 
positive rational number k is the extra input parameter.) 
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Technically speaking, the decision version of MIS is NP-complete, and the optimization version of MIS 
is NP-hard. Since the decision problem and the optimization problem are polynomially transformable to each 
other, in the following, we will simply refer to the MIS problem, without explicitly mentioning which version 
(optimization or decision). In general, which version should be clear from the context. 

One way to solve the MIS problem is by maximizing a quadratic binary function y (also known as pseudo- 
boolean function) defined in the following theorem (Theorem 5.1 in Q3J). 

Theorem 2.1. If Jij > min{cj, Cj}for all ij G E(G), then the maximum value of 

y{x\, . . . , Xn) — ^ C-iXi ^ JijXiXj (2) 
ieV(G) U'6E(G) 

is the total weight of the MIS. In particular if Jij > minjcj, Cj} for all ij G E(G), then mis(G) = {i G V(G) : 
x* = 1}, where (x\, ...,x* n )= argmax (llj ..., Xn ) 6 { ,i}« K^i, • • -,x n ). 

Notice that in this formulation, we only require Jij > min{ci, Cj}, and thus there is freedom in choosing this 
parameter. It is easy to see (by changing the variables Xi = ^4p) that MIS is equivalent to minimizing the Ising 
energy function : 

£(si,...,s n ) = hiSi+ JijSiSj, (3) 

ieV(G) ijeE(G) 

which is the eigenf unction of the following Ising Hamiltonian: 

H C = hia i + J v a i a 3 (4) 

ieV(G) ijeB(G) 

where h{ = Ejenbr(i) J n ~ 2c i' nbr (0 = {j '■ U G E(G)}, for i G V(G). 

In other words, using the Ising Hamiltonian He as the problem Hamiltonian, we have an AQO algorithm for 
solving the MIS problem. At the end of the evolution, the ground state of He, say \x*X2 ■ ■ ■ x* n ), corresponds to 
the maximum-weight independent set, namely mis(G) = {i : x* = 0}. 

2.2 Exact Cover 

Formally, the Exact Cover is as follows: 

Input: A set of m elements, X = {ci, C2, • • . , c m }; a family of n subsets of X, S = {Si, S2, ■ ■ ■ , S n }, 
where Si C X 

Question: Is there a subset I C {1, . . . , n} such that Ui e jSi = X, where Si n Sj ■ = for i ^ j G I? Here 
{Si : i G J} is called an exacf cover of X. 

See Figure [T] for an example. In particular, in this example, each element a G X appears exactly in three 
subsets. This kind of special instance is referred as EC3, which can then be polynomially reducible (<p) to the 
positive l-in-3S AT problem. 

EC3 <p positive l-in-3SAT. Given an instance of EC 3 with an m-element set X and n subsets Si, S n , 
we construct a 3CNF boolean formula ^f(xi, . . . , x n ) = C\ A . . . A C m with n variables and m clauses. For each 
element c, G X, let S^, Si 2 , Si 3 be the three sets that consist of Cj. By definition, exactly one of these three sets 
needs to be in the exact cover of X. Therefore, if we define a binary variable Xk for each set Sk such that x^ = 1 
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Figure 1: X = {ci, c 2 , c 3 , c 4 , c 5 }, and 5 = {Si, S 2 , . . . ,S 7 }, with Si = {ci,c 2 ,c 4 }, S 2 = {ci,c 2 ,c 5 }, S 3 = 
{ci,c 3 ,c 4 },S 4 = {c 2 ,c 3 }, S 5 = {c 3 }, S 6 = {c 4 ,c 5 }, S 7 = {c 5 }. Here {Si,S 5 ,S 7 } is the exact cover of 
X. In this example, each element a appears in exactly three subsets. For instance, element c\ (the red ball) 
appears in Si, S 2 and S 3 . (a) The problem instance can be viewed as a graph Gm with n vertices, where vertex 
i corresponds to the set Si, and there is an edge between two vertices if and only if Si and Sj overlap (share a 
common element). There is an exact cover of X if and only if the weight of mis(Gj\/) is m > where the weight 
of vertex i is the number of elements in Sj. (b) The problem can also be reduced to the following instance of 
positive l-in-3SAT: V(xi, . . . , x-j) = C\ A C 2 A C 3 A C 4 A C 5 where C\ = x\ V x 2 V x 3 , C 2 = x\ V x 2 V x 4 , 
C 3 = x 3 V x 4 V X5, C 4 = »i V x 3 V X6> C5 = x 2 V X6 V X7. There is an exact cover of X if only if is satisfiable 
in that there is exactly one variable in each clause is true. 



if Sk is in the exact cover, and Xk = otherwise, we require exactly one variable in d = x^ V Xi 2 V Xj 3 to 
be satisfied. That is, there is an exact cover to the original problem if and only if the formula ^{xi, . . . , x n ) = 
Ci A . . . A C m is satisfiable in that there is exactly one variable in each clause is satisfied, which is the positive 
l-in-3SAT problem. See Figure[jjb) for an example of the reduction. 

One way of solving the positive l-in-3SAT problem is to minimize the following clause- violation cost func- 
tion 
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= 1 



is satisfiable if and only if the minimum of Sq, is zero (i.e. no violation). The corresponding problem Hamil 



tonian based on this cost function as used by Altshuler et al.[9] (and Young et al. lfl6l ^ 



is 



u A = E & 

i£V(G EC ) ii£E(G EC ) 

where B\ is the number of clauses that contains variable Xi, and ly is the number of clauses that contains both 
Xi and Xj, and V(Gec) = {1, • • • , n}, and E(Gec) = {ij '■ %i and Xj appear in a clause.}. 

Next, we show that there is another way of looking at the Exact Cover problem, which leads to a simple 
polynomial reduction from Exact Cover to MIS. 

Exact Cover <p MIS. Given an instance of Exact Cover with an m-element set X and n subsets Si, ■ ■ ■ , S n . 
Recall that each element in X can appear in exactly in one set. If two sets overlap, e.g. Si and S 2 in Figure [T] 

'The sign of of term is in opposite because they use Xi = instead, iy was called Jy(= |(Jy + Jji)) in H). 
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then they can not both appear in the exact cover. Therefore, if we construct a graph Gm with n vertices, where 
vertex i corresponds to the set Si, and there is an edge between two vertices if Si and Sj overlap. Then an exact 
cover of X will correspond to an independent set of Gm- Moreover, if we let the weight of vertex i be the number 
of elements in Si, then there is an exact cover to the original problem if and only if the weight of mis(GA/) is m. 

For EC3, it is easy to see that Gec and Gm are exactly the same because there is one-one corresponding 
between the variable Xi and the set Si (V(Gec) = V(Gm)), and and Xj appear in a clause" is equivalent 
to "Si and Sj share a common element" (E(Gec) = E(Gm))- Based on this reduction, we therefore have the 
following problem Hamiltonian for the same problem: 

U o = E E Jii-2B i )v!+ E J v<°i < 6 > 
ieV(G EC ) \jenbr(i) J ijeE(G EC ) 

where > min{i?j, Bj}. 

Comparison of Ha and He- For comparison, recall that 2Bi = X^enbr(i) A?> ^ et us wrrte Jij = ^hj + Dij, 
where Dij > mm{Bi, Bj} — 2/jj. Then we have 

h c = m A + E E *Vf+ E D ^ a J- 

ieV(G EC )ienbr(i) ijeE(G EC ) 

Alternatively, as we show in 1171 . for > 0, we can view He as obtained from the following clause- violation 
cost function: 

m 

. • • ,x n ) = ^(x kl + x k2 + x k3 - I) 2 + ^2 DijXiXj (7) 

fc=l ijeE(G EC ) 

= £y(xi, . . . ,x n ) + D^XiXj (8) 

ii6E(G EC ) 

Notice that the extra term (^ijeE(G EC ) Dij x i x j) wn l not contribute to the energy function for the truth as- 
signment because at least one of {xi, Xj} (for ij G E(Gec)) W1 H De zero. In other words, both and are the 
cost functions for the same problem instance ^, where > can be arbitrary. The question is: does this extra 
term (for arbitrarily chosen Dij > 0) matter to the performance of the adiabatic algorithm? 

In [9], Altshuler et al. claimed that the AQO algorithm with problem Hamiltonian Ha failed with high 
probability for randomly generated instances of EC3 due to the Anderson localization(AL). The authors claimed 
that the correctness of their argument did not rely on the specific form of the problem Hamiltonian for EC3, 
but only depended on the properties of the problem instance Bi and Jy. However, our problem Hamiltonian 
Tic belies their claim. Their argument necessarily depends on the energy function of the problem Hamiltonian. 
While the energy function for Ha only depends on Bi and Jy , the energy function for He also depends on the 
extra whose values have a range to choose. Therefore, their argument which depends on the random property 
of parameters does not apply to the algorithm based on He because are not random. Perhaps, their main 
message is that from the physics point of view, if AL occurs, then the AQO algorithm will require exponential 
time, as shown by Amin and Choi in |fT9l . So what are the necessary conditions for AL to occur? Can the 
arguments in Q be modified such that AL still occurs, and the AQO algorithm with He for arbitrarily chosen 

will also fail? Here we emphasize that while the connectivity of the graph is random, the values of D^ (at 
least min{Bi, Bj} — 2Iij) are free to choose. To show that the AQO algorithm with He fails, one would need 
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to show that the algorithms correspond to all possible values of D^ all fail. Or alternatively, can one show that 
there does not exist Dij (not necessarily efficiently computable) such that the corresponding adiabatic algorithm 
is not at least exponential? In |[T6ll . Young et al. used QMC to show that the minimum spectral gap g m - m of the 
AQO algorithm (with the same problem Hamiltonian Ha) is exponentially small. It will be interesting to see the 
g min result (for the same set of instances) using this new problem Hamiltonian He f° r some D^. 



2.3 3SAT 

Similarly, for 3SAT, there is a well-known reduction to MIS (which is one of the first NP-complete reductions, 
to show the NP-hardness of MIS) (H. For completeness, here we recall the reduction: 



3SAT <p MIS. Given a 3SAT instance ^(xi, . . . , x n ) = C\ A . . . A C rn with n variables and m clauses, we 
construct a (unweighted) graph Gsat as follows: 

• For each clause Q = y^ V yi 2 V yi 3 , we construct a triangle with three vertices labeled accordingly, i.e., 
with yi 1 ,yi 2 ,yi 3 , where yj £ {xj,Xj}. Therefore, Gsat consists of 3m vertices. 

• There is an edge between two vertices in different triangles if there labels are in conflict. That is, for i / j, 
i s jt G E(Gsat) if and only if y is = y jt . 

One can then show that ^ is satisfiable if and only if Gsat has a MIS of size m. See Figure [2]for an example. 




Figure 2: Graph Gsat corresponds to . . . , X5) = (xiVx2Vx3)A(x2Vx3Vx5)A(a;iVx4Vx5) A(xiVx3\/x5). 

^ is satisfiable if and only if Gsat has a MIS of size 4. 



In ifTOl . van Dam and Vazirani constructed a special set of 3SAT instances. They showed that the clause- 
violation cost function based AQO algorithm would take exponential time for these instances because of the 
exponential small g mm . Farhi et al.[18] showed that the exponential small gap could be overcome by different 
initial Hamiltonians. Here we point out also that the argument in iPTOl does not apply to the MIS based problem 
Hamiltonian. In order to show that the exponential small gap, their argument requires the cost function to be 
"e-deceptive monotone". For completeness, here we recall the definitions: 

Definition 2.2. (Definition 1 in H10\l ) Consider the hypercube {0, l} n with the partial ordering on its strings 
x X y if and only if Xi < yifor all i = 1, . . . , n. A function f : {0, l} n — > R is monotonically decreasing if 
x H y implies f(x) > f(y). Similarly for monotonically increasing functions. 

Definition 2.3. (Definition 2 in H10\l ) A cost function f : {0, l} n — > M is "e-deceptive monotone" if it is 
monotone for all but a fraction e vertices comprising the top layers of the hypercube. 



However, it is easy to see that the cost function (i.e. the pseudo-boolean function y in Theorem 2. 1 ) for MIS is 



not monotone. For example, 3^(1, 0, 0, . . . , 0) > y(l, 1, 0, . . . , 0) where vertex 1 and vertex 2 are adjacent, while 
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y(l,0, 0, . . . , 0) < 3-7(1, 0, 0, . . . , 1) where vertex 1 and vertex n are independent. Therefore their argument no 
longer applies to this cost function, and one can no longer conclude that this AQO algorithm also fails for the 
instances they constructed. 

3 Discussion 

AQO was originally proposed in j3j H as an energy minimization algorithm that aims to use coherent quantum 
evolution to avoid trapping in the local minima that trip classical algorithms of NP-hard optimization problems. 
Early on, van Dam et al. lITOlfTTI and Reichardt lPT2l showed that the AQO algorithm failed to avoid local minima 
and would take exponential time for computing the minimum of some cost function /, where their AQO algorithm 
is defined such that / is the energy function of the problem Hamiltonian. In other words, in their formulation, the 
energy function of the problem Hamiltonian of their AQO algorithm is same as the cost function of the problem. 
However, given a problem, there are in general many possible cost functions with the same minimum. Thus, 
their AQO algorithm is just a specific algorithm of our more general AQO definition. This perhaps rather obvious 
difference unfortunately is not clear to even quantum computation experts who are not working on the adiabatic 
quantum computation, and the results of van Dam et al. lfTOl ITTTl and Reichardt[ 12] were widely believed that 
they held for all AQO algorithms. Technically speaking, the lower bound of this specific AQO algorithm for 
computing the minimum of / does not constitute a lower bound for all AQO algorithms of the same problem. 
Nevertheless, do these results provide "convincing evidence" that AQO would fail to solve problems with many 
local minima? For this purpose, we constructed a family of MIS graphs in which there are exponentially many 
local minima, our initial results, which were explained by the first order quantum phase transition lTT9ll . agreed 
with the speculation - that is, the system got trapped in the local minima and the particular AQO algorithm 
failed. However, in Ifl4l . we showed that the exponential small gap (caused by the system getting trapped in the 
local minima) could be overcome by changing the parameters of the problem Hamiltonian without changing the 
problem (that is, the global minimum and the exponential local minima remain the same). Although the result 
is only numerical and supported by visualization, this small example serves to clarify that it is not sufficient to 
consider one specific problem Hamiltonian (and thus specific AQO algorithm) for proving the adiabatic lower 
bound of a problem. 

After the appearance of the arXiv version of this paper, F. Krzakala pointed out to the author that the failure 
of "specific" AQO algorithms for the random instances of some NP-complete problems was discussed in their 
work ll20l l2TTl . Note that also, some initial Hamiltonian, such as the projection Hamiltonian used in (22|, can 
make the corresponding AQO algorithm fail as shown in ll23t . Recently, it was also shown in ll24l that the lack of 
structure in the cost function of the problem Hamiltonian can make the corresponding AQO fail. Farhi et al. ll25ll 
suggested that an adiabatic quantum algorithm should be run on each instance with many different random paths 
(as an integral part of AQO). 

In this paper, we describe different AQO algorithms for the NP-complete Exact Cover and 3SAT problems. 
These algorithms are rather straightforward, from the computer science point of view, as they are based on the 
simple NP-complete reductions. However, they serve to further clarify the distinction between one specific prob- 
lem Hamiltonian AQO algorithm and general AQO algorithms for a problem. In particular, we make clear that 
the arguments in both van Dam and Vazirani [10] and Altshuler et al. [9J for their specific AQO algorithm do not 
apply to the AQO algorithms we describe here. While from the physics point of view, "Anderson Localization" 
makes an AQO algorithm fail, does AL necessary occur for all AQO algorithms of the same problem? What are 
the assumptions or conditions for AL to occur? In particular, can one modify the arguments in Altshuler et al. 
such that the AQO algorithm described here also fails? Whether the time complexity of these different AQO 
algorithms is exponential or not requires more rigorous analysis. Our counter-arguments are so simple that one 
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does not need to fully understand the arguments in ||9j [TOl nor to have advanced algorithmic knowledge in order 
to understand. This perhaps highlights the challenge of this interdisciplinary research, and calls for more serious 
investigation. 
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